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— Analysis  of  thunderstorm  turbulence  hazardous  to  aircraft  operation  and 
coordinated  Doppler  radar  observations  indicate  a high  potential  for  Doppler 
radar  utilization  particularly  the  mean  velocity  spectrum  breadth  observations 
in  defining  severe  turbulence  areas.  The  Mean  Velocity  Processor  (MVP,  the 
first  real-time  display  of  Doppler  radar  data)  and  the  Multi-moment  Ling 
Display  (MMD),  both  developed  at  NSSL,  are  utilized  with  the  radars  to  study 
vortex  motion,  turbulence,  and  wind  shear  areas.  In  addition,  the  Plan  Shear 
Indicator  (PSI)  developed  by  the  Air  Force  Cambridge  Research  Laboratory  (AFCRL) 
was  also  employed.  A number  of  severe  convective  storms  were  penetrated  by  an 
instrumented  aircraft  directed  into  areas  which  analysis  inferred  to  be  turbu- 
lent. Aircraft  recorded  turbulence  and  concurrent  Doppler  data  are  compared. 
Utilization  of  the  spectrum  breadth  calculated  from  the  mean  vel^ity  data  as 
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PREFACE 

This  final  report  relates  to  the  investigation  of  Doppler  radar  techniques 
for  real-time  display  of  wind  velocities  in  storms  and  measurements  of  param- 
eters which  correlate  well  with  turbulence  hazardous  to  aircraft  operation. 

This  report  presents  the  state-of-the-art  in  Doppler  radar  as  it  applies  to 
convective  weather  and  aircraft  operations. 

A work  summary  leading  to  a real-time  display  of  Doppler  measured  velocity 
data  is  presented.  Detailed  also  is  the  coordinated  measurements  of  thunder- 
storm turbulence  using  aircraft  and  Doppler  radar.  Mesocyclone  and  tornado 
vortex  recognition  signatures  are  described. 

Meteorological  Doppler  radar  possibly  represents  the  most  significant 
new  capability  for  identification  and  avoidance  of  dangerous  weather  condi- 
tions associated  with  thunderstorms  since  introduction  of  standard  weather 
radar. 


The  work  presented  herein  represents  the  cumulative  effort  of  the  staff 
at  NSSL  and  cooperating  agencies.  Particular  mention  is  made  of  Dr.  Richard 
Doviak,  Mr.  Dale  Si  mans  and  Mr.  Larry  Hennington  for  their  contribution  to 
the  Doppler  and  real-time  displays  and  Messrs.  Donald  Burgess,  Rodger  Brown 
and  Leslie  Lemon  in  the  vortex  identification. 
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APPLICATION  OF  DOPPLER  WEATHER  RADAR  TO 
TURBULENCE  MEASUREMENTS  WHICH  AFFECT  AIRCRAFT 


1 . INTRODUCTION 
1.1  Nature  of  the  Problem 

Thunderstorms  represent  a major  hazard  to  aircraft  operations.  Because 
of  thunderstorms,  large  sections  are  removed  from  air  space  otherwise  avail- 
able to  aircraft,  with  accompanying  disruption  of  service  and  increased  cost 
and  energy  consumption.  At  times,  airports  may  be  closed  to  traffic  for  tens 
of  minutes  or  aircraft  detoured  hundreds  of  miles  to  avoid  thunderstorm 
systems.  FAA  Circular  00-24  addresses  problems  of  thunderstorm  avoidance. 

Use  of  improved  diagnostic  techniques  for  warning  of  severe  thunderstorm 
hazards  would  improve  flight  safety  and  increase  usable  air  space. 

1.2  Objectives  and  Scope  of  Study 

The  development  of  Doppler  radar  offers  the  first  practical  method  for 
measuring  wind  fields  in  detail  inside  storms  and  for  virtually  continuous 
profiling  of  the  horizontal  wind  at  various  altitudes.  Identifiable  wind 
features  should  be  related  to  turbulent  areas  and  other  hazards  such  as  wind 
shear.  This  study  emphasizes  thunderstorms  and  has  objectives:  (1)  to 
develop  methods  for  real-time  display  of  Doppler  wind  data  and  other  parame- 
ters in  a manner  easily  interpreted  and;  (2)  to  determine  hazards  represented 
by  these  data.  Signal  processing  and  display  techniques  were  developed  and 
utilized  with  coordinated  aircraft  measurements  of  suspected  turbulent 
areas.  The  following  sections  relate  progress. 

2.  DOPPLER  RADAR  SIGNAL  PROCESSING 

Radar  echo  intensity  and  Doppler  velocity  must  be  considered  together  in 
analysis  and  application. 

Recent  advances  in  integrated  circuits  and  accompanying  reduction  of 
cost  for  digital  equipment  has  led  to  NSSL's  present  digital  integrator  for 
echo  intensity  (reflectivity)  measurements  (Fig.  1-1).  This  is  generally 
superior  in  performance,  more  stable,  simpler,  more  reliable,  and  less  costly 
than  analogue  systems  now  in  operational  use.  In  addition,  the  digital 
system  offers  more  versatility  in  data  handling,  with  output  data  fed  to  both 
a magnetic  tape  recorder  and  to  a PPI  scope  for  real-time  interpretation 
(Fig.  1-2). 

2.1  Mean  Velocity  Processor 

Real-time  display  of  Doppler  radar  velocity  data  has  also  been  inves- 
tigated. 


To  Display 
Equipment 


To  Digital 
Recorder 


System  Timing 


Figure  1-1.  Digital  integrator  block  diagram. 


Figure  1-2.  Contoured  display  of  the 
WSR-57  radar  echo  (>10  dBz)  of  the 
Davis  Storm  at  1700  CST.  Antenna 
elevation  is  0 deg. 


The  "Mean  Velocity  Processor" 
(MVP)  [originally  called  "Octant 
Change  Counter"  (OCC)]  (Sirmans, 
1973)  provides  an  estimate  of  the 
mean  velocities  of  meteorological 
targets  simultaneously  at  200  range 
locations  along  a radial.  It  is 
the  first  device  in  the  world  to 
provide  a contour-mapped  PPI 
display  of  meteorological  target 
velocities  in  real  time.  The  MVP 
processing  rate  is  commensurate 
with  the  radar  output  data  rate, 
and  mean  velocity  estimates  are 
acquired  in  a time  limited  only  by 
the  desired  estimate  accuracy. 

This  system  estimates  mean  veloci- 
ties at  multiple  pulse  volumes 
along  a radial  while  circumventing 
spectral  calculations  and  large- 
scale  storage  by  acting  on  pairs  of 
complex  time  samples  to  compute  the 
echo  scalar  phase  change  over  a 
pulse  repetition  period.  Each 
calculated  phase  is  compared  with 
the  previous  phase  and  the  quan- 
tized displacements  from  consecutive 
pulses  are  stored.  The  consecutive 
displacements  are  summed  algebrai- 
cally for  a fixed  number  of  pulse 
pairs;  the  expected  net  displacement 


Figure  1-3.  PPI  reflectivity 
(A)  and  Doppler  ieotaohe  (B) 
at  2115  CST.  The  elevation 
angle  is  1.9°,  RM  60,  80, 

100  Ton.  Reflectivity  categories 
are  dim  (<21  dBz),  bright 
(21-31),  black  (31-44),  dim 
(44-57),  and  bright  (>57  dBz). 
Velocity  categories  are  dim 
(<13  m sec-l),  bright  (13-21), 
and  brightest  (>21  m sec~^). 
Positive  radial  velocities  are 
angularly  strobed.  Mesocyclone 
type  signature  at  198° /82  km. 


is  proportional  to  the  mean  velocity  in  the  sampling  volume.  A sample 
simultaneous  display  of  mean  echo  intensity  and  Doppler  velocity  estimate 
in  a PPI  mode  is  shown  in  Figure  1-3. 

Doppler  radar  measurements  in  severe  storms  are  significantly  limited 
by  a relationship  between  the  maximum  unambiguous  range  Rg  and  the  maximum 
unambiguously  resolved  ra '’al  velocity  jV^,!  (i.e.  the  Nyquist  limit).  This 
is  defined  by  the  equations  |V^|  = 2(PRF)  X/8  and  Rg  = C/2  PRF  where  C is 
the  speed  of  light,  X is  the  radar  wavelength,  and  PRF  is  the  pulse  repeti- 
tion frequency.  Thus  for  a given  radar,  increasing  PRF  to  give  a greater 
velocity  interval  results  in  decreased  range.  The  study  by  Doviak  and 
Sirmans  (1973)  suggests  use  of  horizontal  and  vertical  polarization  with 
pulse-pair  processing  to  extend  the  range  of  velocities  measurable  with  a 
specified  PRF.  The  proposal  merits  further  study. 

Further  investigation  of  pulse-pair  processing  techniques  is  also 
required  to  obtain  more  complete  information  from  a singly  polarized  system. 
A system  that  calculates  the  vector  phase  change  rather  than  the  scalar 
change  between  successive  pulses  provides  more  accurate  velocity  estimates 
than  the  MVP  when  velocities  are  near  the  Nyquist  limit,  and  is  described  on 
page  6 below. 

2.2  Multi moment  Display 

Doppler  radar  provides,  in  addition  to  the  precipitation  echo  spectrum 
power  (zeroth  moment)--which  is  proportional  to  reflectivity  (Z)— and  the 
mean  Doppler  velocity  (first  moment),  the  velocity  spread  or  spectrum  broad- 
ness, B,  (second  moment  about  the  mean)  of  precipitation  particles  within 
the  pulse  volume. 
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Integration  of  these  three  spectrum  characteristics  into  a single  display 
provides  superior  real-time  information  for  decision  making. 

The  display  developed  in  1974  is  described  by  Burgess  et  (1976)  as 
follows: 

Display  Description 

A minicomputer-graphic  display  terminal  has  been  interfaced 
to  a National  Severe  Storm  Laboratory  (NSSL)  10  cm  Doppler  radar 
having  characteristics  given  by  Brown  et  (1971).  The  algo- 
rithms necessary  to  estimate  Doppler  moments  and  display  controls 
are  accomplished  through  use  of  the  minicomputer.  Magnetic  tape 
also  records  displayed  data. 

To  simultaneously  present  the  three  principal  Doppler  moments 
for  each  pulse  volume,  a field  of  arrows  is  displayed  where  arrow 
length  is  proportional  to  the  log  of  received  power,  arrow  direc- 
tion to  velocity  and  arrowhead  size  to  Doppler  spectrum  width. 

For  examples  presented  here,  arrowhead  size  is  programmed  to  be 
proportional  to  Doppler  wind  shear,  a significant  contributor  to 
spectrum  width. 

Figure  2-1  illustrates  the  above  relations.  Zero  velocity 
is  a horizontal  arrow  pointing  right  and  nonzero  velocities  are 
proportional  to  the  angle  between  the  arrow  and  its  zero  posi- 
tion. Clockwise  rotation  from  the  zero  position  denotes  negative 
velocities  (toward  the  radar)  and  counterclockwise  rotation 
denotes  positive  velocities  (away  from  the  radar).  The  hori- 
zontal arrow  pointing  left  corresponds  to  the  maximum  unambiguous 
velocity  (±34  m s"')  resolved  by  the  radar.  As  the  velocity 
increases  beyond  ±34  m s-1,  the  arrow  rotates  smoothly  and  appears 
as  a lower  velocity  of  opposite  sign  (e.g.,  38  m s"l  appears 
as  -30  m s"^).  Arrowhead  width  reflects  the  largest  absolute 
velocity  difference  between  the  sample  volume  under  consideration 


Figure  2-1.  Relationships  between 
display  arrows  and  the  three 
Doppler  moments.  Arrow  length 
is  proportional  to  received  power ^ 
arrow  direction  to  velocity  and 
arrowhead  size  to  Doppler  spec- 
trum width  or  velocity  difference. 
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and  the  four  surrounding  volumes.  Due  to  the  limited  unambiguous 
velocity,  the  velocity  difference  is  computed  by 

IVt-V,!  when  |V,-V,1^34  m s"^ 

AV  = ' ^ T 

68-lV^-V2l  when  | V^-V2l^34  m s' ' 

to  eliminate  the  potential  for  large  numbers  of  erroneous  differences 
across  the  maximum  unambiguous  'Velocity.  Therefore,  the  difference 
between  +30  and  -28  m s-1  is  displayed  as  10  m s'l,  not  58  m s-1. 
Unfortunately,  real  velocity  differences  >34  m s"^  are  ignored  in 
arrowhead  size  but  arrow  direction  may  still  be  used  to  correctly 
interpret  large  velocity  differences. 

The  display  (Fig.  2-2)  is  a B-SCAN  (range  vs.  azimuth) 
presentation  of  Doppler  moments  where  ranges  in  kilometers  and 
azimuths  in  degrees  are  given  along  the  left  and  bottom  margins, 
respectively.  Top  headers  show  time,  radar  elevation  angle, 
kilometer  range  interval  (dR)  and  kilometer  azimuth  interval 
(dAZ)  between  whole  degrees.  The  display  sector  is  limited  to  15 
X 16  pulse  volumes  but  range  and  azimuth  spacing  between  dis- 
played pulse  volumes  is  variable  and  the  sector  limits  can  be 
changed  quickly.  Thus  it  is  possible  to  check  large  storm 
regions  for  severe  storm  signatures.  A variable  minimum  power 
threshold  is  used  to  separate  noise  from  regions  where  velocity 
estimates  may  be  obtained.  For  the  examples  shown,  a signal-to- 
noise  ratio  of  3.5  dB  is  used  as  the  boundary  for  displayed  data. 


Figure  2-2  Multimoment  display 
for  Stillwater  storm  of  13  June 
19?5  before  tornado  formation. 
Timet  elevation  angle t range 
interval  (dR)  and  azimuth 
interval  (dAZ)  at  the  nearest 
range  between  arrows  are 
given  across  display  top. 
Mesoayolone  center  is  indi- 
cated by  (+). 


Velocities  are  computed  using  the  pulse  pair  technique  to 
determine  spectral  moments  from  signal  covariance  (Berger  and 
Groginsky,  1973).  To  reduce  computational  time,  alternate  sample 
pairs  (e.g.,  samples  1 and  2,  5 and  6,  etc.)  are  processed  to 
obtain,  from  64  uniformly  spaced  samples,  16  pulse  pairs  to 
estimate  Doppler  moments.  The  standard  error  of  velocity  esti- 
mates is  about  0.7  m S"1  for  spectra  width  equal  to  3 m s-1 , a 
typical  value.  The  use  of  alternate  pairs  results  in  a negli- 
gible increase  of  variance  over  that  predicted  by  Berger  and 
Groginsky  if  the  full  63  sample  pairs  were  processed. 2 These 
errors  are  within  the  accuracy  expected  when  interpreting  the 
velocity  displays  and  are  considered  reasonable  for  real  time 
use. 

This  display  is  very  effective  and,  beginning  in  1975,  is  used  in  the 
decision  making  processes  whenever  penetration  aircraft  are  flying. 

2.3  Plan-Shear  Indicator 

A third  type  of  display,  the  Plan-Shear  Indicator  (PSI)  developed  by 
the  then  Air  Force  Cambridge  Research  Laboratory  (AFCRL)  (Armstrong  and 
Donaldson,  1969)  was  also  tested.  This  is  discussed  in  detail  by  Lee  and 
Kraus  (1975).  The  PSI  display  (Fig.  2-3)  indicates  shear  by  a differential 
displacement  of  adjacent  precipitation  targets  producing  "wiggles"  in  the 
display.  The  width  of  the  concentric  area  is  proportional  to  the  second 
movement.  However  the  scale  factor  did  not  permit  use  of  this  characteristic. 
Results  using  these  displays  are  detailed  in  the  next  section. 

3.  PLAN-SHEAR  INDICATOR  AND  AIRCRAFT  MEASUREMENTS 
OF  THUNDERSTORM  TURBULENCE 

A Doppler  experiment  Spring  1973,  in  Oklahoma,  involved  Air  Force 
Cambridge  Research  Laboratories  (AFCRL,  presently  Air  Force  Geophysical 
Laboratory),  Aeronautical  Systems'  Division  (ASD),  Federal  Aviation  Admin- 
istration (FAA),  Colorado  State  University  (CSU)  and  National  Severe  Storms 
Laboratory  (NSSL). 


3.1  Test  Operation 

The  NSSL  Doppler  radar,  Norman,  Oklahoma  was  modified  to  operate  in 
either  the  PSI  or  standard  mode  since  the  differing  requirements  on  pulse 
repetition  frequency  of  the  two  systems  did  not  permit  simultaneous  operations. 
During  a portion  of  the  program,  the  radar  was  alternated  between  the  two 
modes  for  a series  of  penetrations  by  instrumented  aircraft. 

Ground  clutter  out  to  38  km  and  second  trip  considerations  limited 
alternate  sampling  operations  to  ranges  from  38  km  to  114  km  and  from  152  km 
to  164  km,  substantially  reducing  the  number  of  opportunities  for  comparisons. 


Private  communication,  Dr.  Dusan  Zrnic,  a visiting  scientist  at  the 
National  Severe  Storms  Laboratory. 
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When  a storm  with  intensity  larger 
than  30  dBZ  approached  or  appeared 
within  the  acceptable  ranges, 
aircraft  were  launched  and  vectored 
to  the  storm  area  avoiding  reflecti- 
vities zones  greater  than  50  dBZ 
where  damaging  hail  is  common 
(Burnham  and  Lee,  1969). 

The  PSI  Doppler  scope  was 
viewed  for  indications  of  shear. 

When  a shear  area  was  observed, 
(wiggles  present)  the  aircraft  was 
vectored  through  the  area.  Aircraft 
tracking  was  accomplished  with 
transponder  data  fed  into  the 
WSR-57  10-cm  weather  radar  system 
and  displayed  simultaneously  with 
contours  of  radar  echo  reflectivity. 
This  photographed  display  provides 
an  aircraft  position  every  20  seconds. 
For  analysis,  straight  line  interpo- 
lation was  used  between  each  recorded 
position  (Fig.  3-1).  WWV  time 
signals  were  used  to  coordinate 
aircraft,  weather  radar,  Doppler 
radar,  and  voice  data.  A 1 km 
radius  circle  of  error  is  probable. 
Aircraft  data  were  recorded  at  50 
per  second  and  a five-point  smoothed 
average  provided  0.1  s values  for 
computation  of  true  vertical  veloci- 
ties (w)  and,  as  a measure  of 
turbulence,  derived  gust  velocities 
(Ude)  (Hombolt,  et  ai,  1964).  The 
penetration  true  airspeed  was 
approximately  170  m s"‘  with  the 
F-lOO,  210  m s-1  with  the  F-101. 

Thus  0.1  s data  points  correspond 
to  observations  about  every  20  m. 

The  PSI  has  a range  resolution  of 
855  m and  the  0.8°  beam  width 
provides  an  azimuthal  resolution  of 
about  1.1  km  at  80  km.  Hence,  we 
are  concerned  with  turbulent  encoun- 
ters of  several  seconds  or  more. 

More  specifically,  the  0.1  s aircraft 
data  output  corresponds  to  a maximum 
unambiguous  frequency  (f^a^) 

5 Hz  and  a corresponding  minimum 
) along  the  flight 
by 


wavelength  (1^,^ 
path  represented 


target  moving  toward 


RADAR 


Figure  2-3.  PSI  display  for 
stationary  targets  (left)  and  a 
moving  target  (right).  The  moving 
target  is  located  at  the  same  dis- 
tance from  the  radar  as  the  nearest 
stationary  target^  but  is  displaced 
from  it  on  the  PSI  display  by  an 
increment  which  depends  on  its 
velocity  (Armstrong  and  Donaldson, 


1969). 


Figure  3-1.  16  June  1973,  WSR-57 

weather  radar  reflectivity  iso- 
echo contours  display  with 
aircraft  transponder  beacons 
superimposed.  Point  A is  the 
beacon  return  from  the  F-lOO 
at  1357:25  CST;  the  dotted  line 
indicates  aircraft  path.  Range 
marks,  at  40  km  intervals. 


where  V is  aircraft  speed  and  V the  wind  speed  along  the  flight  path. 
Assuming  the  winds  are  contained  within  the  maximum  unambiguous  velocity  of 
the  Doppler  radar  (±34  m s"l)  and  that  the  Taylor  hypothesis  holds,  L^in  is 
30  to  40  m.  The  maximum  wavelength  resolvable  is  equipment  dependent  and  is 
estimated  to  be  about  2000  m (Ryan  et  al^. , 1971). 

In  the  Doppler  data  field,  Lmin  is  dependent  on  range  and  aspect  angle 
due  to  radar  beam  divergence  and  differences  between  range  and  angular  reso- 
lution. In  the  series  of  complete  Doppler  spectrum  data,  the  range  sampling 
interval  of  300  m yields  L^in  = 600  m if  the  aircraft  flies  along  a radar 
beam.  At  the  working  distance  of  60  km,  the  0.8°  beam  width  provides  Lmin 
of  about  1600  m.  Data  obtained  within  a pulse  volume  are  the  product  of 
several  integrated  factors  throughout  the  total  volume  which  is  quite  large 
(3.308  X 108  m3  at  60  km)  and  it  is  difficult  to  resolve  turbulent  wave- 
lengths within  the  volume.  Thereby  we  see  that  pattern  correspondence 
rather  than  point-to-point  correspondence  must  be  utilized  in  correlation 
studies. 

Another  consideration  is  that  the  Doppler  primarily  measures  the  hori- 
zontal wind  component  at  low  elevation  angles  used  for  PPI-type  displays. 
With  aircraft,  the  vertical  component  is  measured.  However,  Houbolt  et  al . , 
(1964),  Ashburn  ^al.,  (1970)  and  others  have  shown  that  turbulence  affect- 
ing aircraft  above  tRie  boundary  layer  is  mainly  isotropic,  particularly  so 
in  convective  clouds.  This  we  assumed  in  the  present  study. 

The  F-101  made  penetration  on  6 days  (Table  3-1);  not  all  data  from 
these  flights  are  available.  The  F-lOO  aircraft  penetrated  thunderstorms  on 
five  days.  On  three  of  these  days  (12  penetrations)  it  was  possible  to 
compare  PSI  presentation  and  Ude*  Derived  gust  velocities  were  divided  into 
four  categories  as  shown  in  Table  3-2.  When  the  aircraft  recorded  only 
light  turbulence  or  less,  the  PSI  arcs  were  smooth,  for  example  15  June 
F IRl  and  16  June  F 1R5  (Table  3-3). 


TABLE  3-1 

Aircraft  Operation  Days  when  PSI  Data  Recorded 


F-101 

F-lOO 

26  April  1973 

1 June  1973 

6 May 

4 

21 

11 

23 

15 

26 

16 

1 June 

* ' > • - 


TABLE  3- 

•2 

Interval 

Approximate  Turbulence  Class  (NASA) 

3.0  < Ujjg  < 6.0  m s'^  (10-19  ft  s'^ 

Light 

6.1  i Ujg  s 9.1  m (20-29  ft  s’b 

Moderate 

9.2  < Ujg  < 12.1  m s"^  (30-39  ft  s"^ 

Severe 

^ 12.2  m s"^  (^  40  ft  s"^ 

Extreme 

TABLE  3-3 

Date  (1973) 

Time(s) 
in  Turbulence 

and 

Cat.  BlCat.  C 

Flight 

(MDT)  (SVR) 

PSI  Display 

4 June  F 1R5 


15  June  F IRl 


F 2R1 


F 2R3 


16  June  F 1R3 


F 1R4 


F 1R5 


F 1R6 


16  June  F 2R1 


F 2R2 


F 2R3 


F 2R4 


Small  wiqqles 


Wlqqles 


Smooth  (Dim  Photo 
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Other  turbulence  categories  did  not  exhibit  as  good  a correlation. 

Moderate  turbulence  was  encountered  on  ten  penetrations.  Seven  of  these 
also  contained  0.3  s or  longer  embedded  severe  turbulence  areas.  Example: 
on  16  June  the  PSI  indicated  radial  wind  shear  during  penetrations  1 (Fig.  3-2) 
and  3 of  the  second  flight.  In  both  cases,  severe  turbulence  (Fig.  3-3) 
was  encountered  in  the  vicinity  of  the  wiggles.  However,  there  was  no 
indication  of  shear  along  the  flight  path  in  the  moderate  turbulence  area 
preceding  the  severe  turbulence.  For  the  second  run  of  the  same  series 
(F  2R2)  the  PSI  data  appeared  to  show  no  high  shears  in  the  severe  turbu- 
lence areas.  In  the  moderate  turbulence  encounter  areas  there  was  no  evi- 
dence of  wiggles. 

15  June  F 2R1  and  F 2R3  had  areas  in  which  the  One's  indicated  moderate 
turbulence  and  in  F 2R1  a short  burst  of  severe  turbulence.  Along  the  air- 
craft track  the  PSI's  for  both  flights  were  smooth.  A small  area  of  weak 
shear  was  indicated  on  the  PSI  about  10  km  SW  of  the  point  where  the  aircraft 
experienced  some  moderate  turbulence— not  apparently  related  to  the  aircraft 
encounter. 

16  June,  F 1R3  was  associated  with  a smooth  PSI  although  relatively 
long  periods  of  moderate  turbulence,  and  a short  period  of  severe,  were 
recorded.  On  F 1R4  once  again  the  PSI  showed  no  wiggles  along  the  flight 
track  but  did  indicate  some  shear  10  km  NNE  of  the  severe  turbulence  encounter. 

On  F 1R6  of  the  same  day  there  is  an  area  of  shear  along  the  flight 
path  very  close  to  the  location  where  the  aircraft  measured  severe  turbu- 
lence. However,  at  the  time  of  the  PSI  photo  (2  minutes  before  the  turbu- 
lence encounter  time)  the  location  of  shear  is  several  thousand  feet  higher 
than  the  aircraft  location.  Closer  to  the  aircraft  altitude,  but  a minute 
or  two  earlier  indication  of  shear  was  debatable. 

Table  3-2  summarizes  PSI  indicators  and  turbulence  experienced  during 
penetrations  for  all  full  data  set  cases.  An  example  of  vertical  velocities 
encountered  is  shown  by  the  time  history  (Fig.  3-4)  calculated  from  F-101 
data  for  26  May  1973. 


Figure  3-2.  PSI  for  1710  CST  16  June  Figure  3-3.  Penetration  track  time 
1973.  4°  tilt.  20  km  range  marks.  and  turbulence  intensities  encoun- 

tered corresponds  to  flight  in 
Figure  4. 


3.2  Summary 


Moderate  or  severe  turbulence  was  encountered  in  all  cases  when  the  PS I 
displayed  wiggles  along  the  aircraft  flight  path,  but  wiggles  were  not 
present  with  all  turbulence  encounters.  Thus,  it  apoears  from  these  cases 
that  turbulence  up  through  moderate  (Ujjg  < 9.1  m s"')  may  escape  detection 
by  the  PSI.  Where  severe  turbulence  (Ufje  ^ 9-1  m s-l)  repeatedly  was  encoun- 
tered, the  PSI  showed  transient  (less  than  one  minute  duration)  shear  areas 
along  the  flight  path.  Arc  deformations  apparently  have  an  operational 
detectability  threshold  associated  with  wind  shears  > 1.5  x 10"2  s"'. 

4.  DUAL  DOPPLER  AND  SINGLE  DOPPLER  OBSERVATIONS  AND 
CONCURRENT  AIRCRAFT  THUNDERSTORM  PENETRATIONS 

4.1  General 

In  1974  an  F-4-C  aircraft  replaced  the  F-lOO,  a second  Doppler  radar 
was  installed  at  Cimarron  Airport,  Oklahoma  City,  Oklahoma  and  as  noted  in 
Section  2.2  the  multi-moment  display  was  available.  The  second  Doppler  used 
in  conjunction  with  the  Norman  unit  makes  possible  the  calculation  of  winds 
within  a thunderstorm.  The  second  Doppler  also  permits  comparison  of  thun- 
derstorm parameters  as  seen  from  different  viewing  angles.  Operational 
procedures  remained  essentially  the  same  as  detailed  in  Section  3.1.  Average 
penetration  airspeed  for  the  F-4-C  is  150  m s-1.  The  F-4-C  was  on  station 
at  Tinker  AFB,  Oklahoma  (TIK)  from  15  April  through  29  April  1974  and  15  May 
through  30  May  1974.  Eight  flights  were  made  during  this  time;  dates  and 
penetration  altitude  are  shown  in  Table  4-1.  Penetration  altitudes  are  as 
low  as  possible  to  reduce  contamination  of  the  Doppler  velfcity  by  vertical 
velocities  associated  with  falling  precipitation.  In  197t  the  aircraft  was 
available  18  May  through  10  June  and  nine  flights  on  seven  days  were  flown. 
For  1976  the  F-4-C  was  at  TIK  14  May  through  17  June.  Thunderstorm  activity 
was  below  normal  and  only  six  flights,  two  of  which  were  gust  front  investi- 
gations at  4-5,000  ft  MSL,  were  flown.  During  these  days,  45  of  the  penetra- 
tions had  data  sets  which  were  sufficiently  complete  to  be  analyzed. 


FLIGHT  FIOI  PENT  5 MAY  26.  1975 


Figure  3~4,  Vertical  wind  veloci- 
ties during  the  thunderstorm 
penetration  26  May  1973. 

(CEX  = cloud  exit.) 


1458  37  1458  47  1458  57  1459  07  1459  17  1459  27 
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TABLE  4-1 


F-4-C  AIRCRAFT  PENETRATION  FLIGHTS 


1974 

1975 

lATE 

ALTITUDE 

DATE 

ALTITUDE 

•il  21 

12  and  15,000  ft 

May  22 

10,000  ft 

‘il  28 

15,000 

May  24(2) 

8,000 

21 

16,000 

May  26 

12,000 

23 

16,000 

May  29 

15,000 

24 

14,000 

June  8(2) 

20,000 

25(2) 

14,000 

June  5 

(no  data) 

31 

13  and  15,000 



May  19 

(no  data) 

ALTITUDE 


4 10,000  ft 

1 16,000 

9 4 and  5,000 

0 5,000 

12  20,000 

17  15,000 


4.2  Data  Analysis 

The  following  Doppler  radar  derived  parameters  were  considered: 
1.  Radar  reflectivity  {Z){dBZ). 


2.  Reflectivity  gradients  ^ . 

3.  Mean  radial  velocity  (v)(m  s“^). 

4.  Maximum  velocity  gradient  (shear)  x = s \ 

5.  Spectrum  breadth  of  the  mean  velocity  B = m s 

. 9^f  _ .-1  -1 


6.  Gradient  of  the  gradient  = s”  m”  . 

2 -1  -1 

7.  Laplacian  V f = s m . 

when  V = radial  velocity 

r = distance  along  radial 
X = general  distance 
0 = angular  measurement 
c = constant. 

In  the  initial  approach  to  determining  potential  significant  parameters, 
a regression  equation  computer  program  was  developed  in  the  form 

f(x)  = A + A^Y^  + AgYg  ...  + A^Y. 


where  x = derived  gust  velocity  (U^e);  Y],  Y2  ...  Y,  are  independent  variables 
such  as  reflectivity,  mean  wind,  wind  shear  and  A],  A2  ...  are  weighted 
coefficients  related  to  the  correlation.  Using  the  first  two  years  of  data, 
as  input  no  single  significant  element  was  apparent.  We  thus  considered  the 
problem  from  a second  viewpoint:  if  a parameter  is  to  be  a reliable  indi- 
cator, then  it  must  exceed  a threshold  value  whenever  turbulence  above  an 
assigned  value  is  encountered.  This  led  to  a second  computer  program  develop- 
ment which  took  the  form:  if  Uje  > M then  Yi(niax)  - ^i-  That  is,  whenever 
the  derived  gust  velocity  exceeds  a critical  value  (M)  then  the  indicator 
parameter  (Yi)  equals  or  exceeds  a set  value 

This  is  depicted  in  time (space)  cross-sections  of  the  variables  along 
the  flight  path.  The  details  of  the  Doppler  data  processing  program  are 
given  in  Appendix  A and  the  time- hi  story  concept  in  Appendix  B. 

4.3  Relation  of  Doppler  Spectrum  Broadness  to  Aircraft  Turbulence 

2 Figure  4-1  is  a typical  plot.  Values  of  Uje,  B,  3v/3x,  3v/3x3y,  and 
V'^v  are  plotted  using  the  corresponding  linear  scales  at  the  left.  The 


A B C D • 


Fiaure  4-1.  8 June  1975  penetration  4:  Time  (apace)  oroaa-aeation  for  maxi- 

mm  values  recorded  for  each  five  seconds  of  flight  or  corresDondiiw  Doppler 
radtxr  volume  during  penetration.  Derived  gust  velocity  — (f'^ 

trum  breadth  = A(m  e-l)j;  velocity  gradient  = B (s-^  x 1000);  gradient  of 
the  gradient  = C(a-1  rrr^  x 1000).  Solid  line  connects  values  of  spectrum 
breadth  and  dashed  line^  the  derived  gust  velocities.  D = Lapladan;  2 = 
tuo  data  points  at  the  same  place:  3 = three  data  points  at  the  seme  place. 

$ is  used  when  one  or  more  •yariables"  are  plotted  at  the  same  point  as  an 
observed  derived  gust  velocity. 
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scales  were  placed  so  that  what  was  expected  to  be  critical  values  were  on 
the  same  horizontal  line.  Thus  we  have  a B = 4 m s-1  on  the  same  line  as 
tide  = 6.1  m s”1  (20  ft  s"^)  corresponding  to  moderate  turbulence. 

Note  in  this  same  figure  how  well  the  turbulence  trend  matches  the 
trend  in  the  width  plot.  Figure  4-2  and  4-3  are  additional  typical  examples. 
Similar  graphs  were  produced  for  all  45  cases.  It  is  apparent  if  a U^e  ^ 

6.1  m s"l  (20  ft  s"l)  (moderate  turbulence)  value  is  considered  as  critical, 
that  the  spectrum  broadness  successfully  flags  these  occurrences.  The  sta- 
tistic probabilty  of  detection  (POD)  is  described  by  Donaldson  et  al. 

(1975)  and  is  given  by 


POD  = X (X  + Y) 


where  "X"  is  a successful  "forecast"--and  "Y"  is  a miss.  We  apply  this  test 
to  a turbulence  "forecast"  using  the  number  of  encounters  when  Uje  ^ 6.1  m s-1 
and  B > 4 m s-1  as  X and  number  of  occasions  Uje  > 6.1  m s-1  and  B < 4 m s"l 
as  Y. 

A B C D ♦ 

10  30  45  100  50; — » ; 


8 24  36  80  40  I 


6 18  27  60  30 


12  18  40  20 


M 1 


29  MAY  1975 
FLIGHT  7 
PENETRATION  4 
(070475) 


0 0 0 0 0 


40  45 


NUMBER  OF  FIVE  SECOND  INTERVALS 


radar  volume  during  penetration.  Derived  gust  velocity  = *(jt  8~^J;  spec 
trum  breadth  = A(m  8~^)i  velocity  gradient  = B (8~^  x 1000) i gradient  of 
the  aradient  = C(8-^  m~^  x 1000).  Solid  line  connects  values  of  spectrum 


L 


breadth  and  dashed  line,  the  derived  gust  velocities . D = Laplaciani  2 = 
two  data  points  at  the  same  place:  3 = three  data  points  at  the  seme  place. 
$ is  used  when  one  or  more  "variables"  are  plotted  at  the  same  point  as  an 
observed  derived  gust  velocity. 
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Table  4-2  is  a tabulation  of  the  number  of  5 second  intervals  during 
which  Ujjo  > 6.1  m s-1  (20  ft  S”l).  The  76  occurrences  during  45  penetrations 
had  B values  > 4.0  within  1 km  on  71  of  these  encounters,  this  yields  a 
POD  = 0.93.  One  of  the  missed  cases  had  B > 4.0  m s-1  within  2 km  which 

would  raise  the  POD  to  0.95.  In  analyzing  two  additional  cases  it  appears 

from  the  recorded  aircraft  elevator  deflection  that  the  U(jg  values  may  be 
contaminated  by  pilot  input  but  the  amount  of  contamination  has  not  been 
determined. 

Doviak  et  (1976)  have  calculated  the  cumulative  probability  of 
spectrum  broadness  within  the  entire  storm  for  two  tornadic  storms  (Fig.  4-4). 

Note  that  for  B > 4 m s-1  this  probability  is  only  about  30%.  For  non- 

severe  storms  the  probability  should  be  even  less. 

It  is  recognized  that  the  second  moment  reflects  not  only  turbulence 
but  also  wind  shear  and  beam  broadening.  These  have  been  discussed  by 


A B C D » 


Figure  4-3.  29  May  1976  penetration  5:  Time  (ejpaoe)  aroea-seation  for  maxi-  I 

rrwm  vofuea  recorded  for  each  five  seconds  to  flight  or  corresponding  Doppler  i 

radar  volume  dxcring  penetration.  Derived  gust  velocity  = *(ft  a-\) ; spec-  ,• 

tmm  breadth  = A(m  velocity  gradient  = B (a-1  x 1000);  gradient  of  j 

the  gradient  = C(a~^  irr^  x 1000).  Solid  line  connects  valves  of  spectrum  { 

breadth  and  dashed  line^  the  derived  gust  velocities.  D = Laplacian;  2 = f 

two  data  points  at  the  same  place:  3 ® three  data  points  at  the  same  place.  j; 

$ is  used  when  one  or  more  'yariables"  are  plotted  at  the  same  point  as  an  ' 

observed  derived  gust  velocity.  I 


TABLE  4-2 


Occurrences  of  all  Derived  Gust  Velocities 

and  Concurrent  Spectrum  Broadness  Values. 


6.1  m s‘^  (20  ft  s"b 


ght 

^nt. 


3-1-75 


3-2-75 

3-3-75 


3-4-75 


3-5-75 


7-2-75 


X = B > 4 (during  interval); 

0 = B < 4 m s'^  N.P.  = no  U 


X = B ^ 4 within  1 km; 
^ 6.1  m s"^  during  run. 


B when 


U^e  ^ 6.1  ms 


-1 


Flight 

Ident. 


de 
B when 


Ujjg  > 6.1  ms 


-1 


X 

7-3-75 

X 

7-1-74 

X 

X 

X 

X 

X 

X 

7-5-74 

X 

X 

X 

8-1-74 

0 

7-4-75 

X 

X 

X 

X 

8-2-74 

X 

7-5-75 

X 

X 

X 

X 

X 

8-5-74 

X 

7-6-75 

X 

X 

X 

5-1-75 

X 

X 

X 

8-6-74 

X 

X 

X 

X 

9-5-74 

X 

5-4-75 

X 

X 

X 

9-9-74 

X 

X 

11-4-75 

X 

10-2-74 

0 (2  km) 

X 

0 (Pilot  Input) 

X 

10-3-74 

0 (Pilot  Input) 

X 

10-10-74 

X 

6-5-74 

X 

X 

6-3-75 

X 

6-4-74 

X 

X 

X 

X 

X 

Flight 

Ident. 


B when 
> 6.1  m s 


Penetrations  which  did  not  Encounter  6.1  ms' 


3-2-74 

9-6-74 

10-11-74 

10-2-75 

6-2-74 

10-4-74 

3-4-75 

10-3-75 

6-6-74 

10-5-74 

4-4-75 

10-4-75 

7-6-74 

10-6-74 

4-5-75 

11-3-75 

9-1-74 

10-8-74 

5-6-75 

Si  mans  and  Doviak  (1973)  in  detail.  This  does  mean  that  there  will  be 
areas  of  high  B which  will  not  be  turbulent.  However  since  these  areas 
comprise  only  a small  portion  of  a storm  an  "overforecast"  appears  to  be 
acceptable. 

The  other  parameters  investigated  were  not  as  successful  and  when 
Quantities  involving  shear  alone  are  used,  the  dependence  on  viewing  angle 
(tangential  shear  cannot  be  determined)  would  appear  to  reduce  the  opera- 
tional effectiveness.  Turbulence  due  to  convective  processes  appears  to  be 
isotropic  (Houbolt,  £t  al.,  1964)  (Ashburn  et  al . , 1970)  and  therefore  the 
spectrum  width  is  less  dependent  on  viewing  an^e. 

We  tried  to  address  the  question  of  high  standard  deviation  locations 
within  a storm.  Insufficient  Doppler  data  are  available  for  any  conclusion 
to  be  final.  However  Figure  4-5  and  4-6  are  presented  to  show  probable 
correlations  between  turbulence  and  the  up-  (down)  drafts.  These  figures 
are  derived  from  a study  by  Ray  (1975)  of  tornadic  storms.  He  uses  the 
equation  of  continuity  and  dual-Doppler  winds  to  calculate  the  3-dimensional 
storm  structure.  We  show  here  only  a 5 km  height  section  for  which  we 
analyzed  the  data.  Note  that  the  reflectivity  maximum  is  north  of  the 
updrafts  (positive).  It  is  also  apparent  that  the  areas  of  large  spectrum 
broadness  are  on  the  edges  of  the  updraft  with  a preference  for  higher 
values  when  a downdraft  is  in  close  proximity.  This  suggests  that  the 
turbulence  is  produced  by  the  shear  (in  the  horizontal)  of  the  vertical  wind 
(Aw/ Ax  or  Aw/ As)  and  that  as  more  advances  are  made  in  thunderstorm  under- 
standing our  ability  to  locate  these  areas  will  be  improved.  In  addition 
vortex  motion  may  also  be  turbulent. 

4.4  Summary 

The  spectrum  broadness  of  the  radial  velocity  as  observed  by  Doppler 
radar  appears  to  be  a useful  indicator  of  thunderstorm  turbulence.  A PDD 
= D.95  is  obtained  for  B > 4.D  m s-1.  These  spectrum  width  estimates  were 
equaled  or  exceeded  in  only  slightly  more  than  3D%  of  the  volume  within  two 


Figure  4-4.  Cumulative  prob- 
ability of  unbiased  spea- 
trum  widthSj  for  echoes 

from  tuio  tomaata  storms. 
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Burgess  (1976)  and  Brown  and  Lemon  (1976)  in  work  partially  supported 
under  this  contract  developed  mesocyclone  and  tornado  vortex  signature 
recognition  criteria  applicable  to  single  Doppler  radar  observations.  As 
mentioned  in  Section  4,  a single  Doppler  radar  measures  only  the  radial 
component  of  motion  and  is  "blind"  to  motions  perpendicular  to  the  beam. 
Thus,  unique  "signatures"  are  desirable  in  identifying  vortex  motion  when 
dual-Doppler  systems  are  not  practical  or  available. 
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Ftgure  4-6.  8 June  19?4:  Speatrum  breadth  (m  s~  ) contours  are  solid 

lines;  vertical  motion  (m  s~^)  contours  are  dashed  lines.  is 

upuard.  Cross  hatched  area  is  downdraft. 
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5.2  Mesocyclone  Signature 


The  mesocyclone  circulation  identification  has  been  summarized  by 
Burgess  (1976)  as  follows: 

If  a single  Doppler  radar  were  to  scan  a stationary  non- 
divergent  mesocyclone  circulation  having  the  tangential  velocity 
distribution  (two  flow  regimes)  of  a Rankine  combined  vortex,  a 
characteristic  velocity  pattern  would  result  (Fig.  5-1).  The 
inner  flow  regime  is  the  vortex  core  (solid  rotation)  with  a 
velocity  distribution: 

C,  where  C,  = 


r = radius 

R = radius  of  maximum  wind 
V = tangential  velocity 

Outside  the  core,  the  second  flow  regime  is  that  of  a potential 
vortex: 


V-r  = where  = V^^-R 


The  mesocyclone  solid  rotation  core  is  easily  detected  and  its 
size  tneasured  in  single  Doppler  data  by  the  distance  between 
velocity  maxima  of  opposite  sign  (see  Fig.  5-1).  The  average 


Figure  5-1.  Single  Doppler  hori- 
zontal mesocyclone  signature  of 
a stationary  Rankine  conjoined 
vortex  (upper)  and  velocity  pro- 
file along  axis  XY  (lower). 
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core  diameter  of  NSSL  data  is  5.7  km  with  diameters  being  larger 
during  early  mesocyclone  life  and  shrinking  during  tornado  pro- 
duction. The  outer  potential  vortex  region  is  difficult  to 

detect  in  single  Doppler  data  but  NSSL  dual-Doppler  data  suggest  |! 

its  circulation  region  is  about  twice  the  core  size.  All  measure-  j] 

ments  for  the  mesocyclone  used  in  this  paper  pertain  only  to 
core  diameter  and  not  to  potential  vortex  diameter. 

The  characteristic  velocity  pattern  or  couplet  results 
because  tangential  flow  perpendicular  to  the  radar  beam  produces 
zero  velocity  along  a radial  line  from  the  radar  through  the 
center  of  circulation.  Closed  velocity  contours  (isodops)  appear 
at  the  radius  of  maximum  wind  where  the  radar  beam  is  parallel  to 
the  tangential  flow.  The  presence  of  convergence  (divergence) 
rotates  or  skews  the  pattern  clockwise  (counterclockwise)  rela- 
tive to  the  radial  line  from  the  radar  to  the  circulation  center. 

Single  Doppler  mesocyclone  signatures  in  tornadic  storms  (Fig. 

5-2)  conform  well  with  the  idealized  model  and  at  most  times  are 
recognizable  at  first  glance.  Signature  vertical  continuity  is 
achieved  throughout  the  storm's  lowest  8 to  10  kilometers.  The 
signature  has  been  verified  as  a closed  cyclonic  circulation 
several  times,  beginning  with  Brown  et  al^.  (1975).  For  the  same 
storm,  Ray  ^ a^.  (1975)  calculated  tRe  mesocyclone  signature 
region  to  be  coincident  with  high-speed  updraft. 

5.3  Mesocyclone  Statistics 

Table  5-1  is  a list  of  the  37  mesocyclones  which  have  been 
confirmed  by  satisfying  the  criteria  previously  mentioned.  Three 
signatures  from  the  original  data  set  had  to  be  eliminated  because 
the  data  collection  mode  failed  to  establish  their  persistence. 

The  37  mesocyclonic  storms  are  a subset  of  the  approximately  350 
storms  scanned  by  single  Doppler  from  Spring  1971  to  Spring  1975. 

At  least  23  mesocyclones  (62%)  are  associated  with  reported 
tornadoes.  Storm  Data  (1971-1975)  was  used  as  the  tornado  report 


•ao*  170*  160*  (50*  19  APRIL  1972 


Figure  5-2.  Single  Doppler  velo- 
city field  with  mesocyclone 
signature  at  the  time  of  a 
tornado  (black  dot). 
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source  (except  for  a few  NSSL  surveys)  and  a measurable  percentage 
of  tornadoes  is  not  contained  in  storm  data,  particularly  those 
occurring  in  rural  areas  after  dark.  Of  the  remaining  14  signa- 
tures, 12  are  associated  with  damaging  wind  or  hail.  Fully  95% 
of  all  mesocyclones  have  produced  some  type  of  surface  damage. 

At  no  time  during  data  collection  did  a verified  tornado  occur 
unless  preceded  by  a mesocyclone  signature.  The  average  lead 
time  before  tornado  occurrence  was  36  minutes.  Some  mesocyclones 
were  well  formed  when  data  collection  began  and  lead  times  could 
have  been  even  longer. 

The  average  mesocyclone  has  a rotational  velocity  of  22  m s"^, 
a horizontal  diameter  of  5.7  km,  and  a vertical  extent  of  7.8  km. 
Very  little  difference  in  rotational  velocity,  horizontal  diameter, 
or  vertical  extent  is  seen  between  tornadic  and  non-tornadic 
mesocyclones.  However,  mesocyclones  which  produce  maxi -tornadoes 
have  smaller  horizontal  diameters,  are  taller  and  rotate  consider- 
ably faster. 


5.4  Tornado  Vortex  Signatures 

The  tornado  vortex  is  especially  hazardous  to  aircraft  and  is  much 
smaller  than  the  mesoscale  system  of  which  it  is  a part.  This  poses  addi- 
tional recognition  problems.  Work  on  this  phase  has  been  reported  by  Brown 
and  Lemon  (1976)  and  excerpts  follow: 

When  one  considers  the  more  typical  (lower  PRF)  Doppler 
radar  with  limited  velocity  resolution,  tornadic  wind  speeds  will 
fold  into  the  unambiguous  velocity  interval.  The  resulting  broad 
spectrum  fills  the  entire  unambiguous  interval,  making  it  extremely 
difficult  to  unscramble  the  peak  tornadic  speeds.  Using  the 
5.4  cm  Air  Force  Cambridge  Research  Laboratories'  Doppler  radar, 

Kraus  (1973)  found  filled  spectra  ( 12.5  m s-1)  in  the  vicinity 
of  the  Brookline,  Massachusetts  tornado  of  9 August  1972.  Zrnic 
et  (1976)  have  estimated  peak  speeds  by  comparing  numerical 
simulations  with  folded  spectra  measured  by  the  10  cm  NSSL  Doppler 
radar. 

For  low  PRF  data  collection  throughout  a convective  storm, 
the  mean  of  the  Doppler  velocity  spectrum  is  the  parameter  of 
primary  interest.  Up  to  the  time  of  the  NSSL  Doppler  velocity 
measurements  in  the  Union  City  tornadic  storm  of  24  May  1973,  the 
mean  Doppler  velocity  field  had  not  been  expected  to  indicate  the 
presence  of  a tornado.  However,  the  data  collected  on  that  day 
did  reveal  a distinct  tornadic  vortex  signature  (TVS). 

The  single  Doppler  velocity  signature  of  a vortex  is  a 
function  of  vortex  size  relative  to  radar  beam  size.  Mesoscale 
vortices— having  dimensions  considerably  greater  than  the  beamwidth— 
produce  a unique  signature  when  mean  Doppler  velocity  values  are 
plotted  as  a function  of  range  and  azimuth  (e.g.,  Donaldson, 

1970;  Burgess,  1976). 
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For  vortices  smaller  than,  but  centered  in,  the  radar  beam, 
previously  cited  literature  suggests  that  the  single  Doppler 
velocity  signature  of  a tornado  vortex  should  be  a broad  velocity 
spectrum  with  a mean  of  zero.  However,  when  the  radar  beam  is 
not  centered  on  the  vortex,  it  is  difficult  to  accurately  predict 
the  spectra  and  mean  Doppler  velocity  fields.  Fortunately,  Zrnic 
recently  developed  a model  that  simulates  a Doppler  radar  looking 
at  a Rankine  combined  vortex  (see  Zrnic  and  Doviak,  1975).  The 
radar  and  vortex  characteristics  are  completely  flexible  such 
that  the  model  can  be  used  to  simulate  Doppler  velocity  measure- 
ments in  vortices  ranging  from  very  small  tornadoes  to  large 
mesocyclones.  The  model  also  allows  the  reflectivity  profile 
across  the  vortex  to  be  varied. 

The  Zrnic  model  aids  understanding  of  data  from  a finite 
Doppler  radar  beam  scanning  across  a vortex.  In  a Rankine  com- 
bined vortex  (heavy  curve  in  Fig.  5-3),  the  tangential  velocity 
increases  linearly  until  the  maximum  velocity  at  the  outer  edge 
of  the  "core"  is  reached  then  decreases  inversely  proportional  to 
radius.  Various  sized  vortices  generated  by  the  model  can  be 
normalized  by  dividing  all  velocities  by  the  maximum  core  velocity 
and  by  dividing  all  lengths  by  the  core  radius.  A uniform  reflec- 
tivity profile  is  used  for  the  simulations  presented  in  this 
paper. 


Figure  5-3,  Theoretical  change  of  mean  Doppler  velocity  azimuthal 
profile  through  a Rankine  combined  vortex  (heavy  line  labeled  0.0) 

08  the  beanuoidth  becomes  progressively  larger  relative  to  the  vortex 
core  radius.  Velocities  and  distances  are  normalized  relative  to 
peak  vortex  tangential  velocity  and  core  radius  of  peak  velocity^ 
respectively. 
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As  indicated  in  Figure  5-3,  when  the  radar  half-power  beam- 
width  is  a small  fraction  of  the  vortex  core  radius  (beamwidth/ 
core  radius  ratios  much  less  than  one),  the  Doppler  velocity 
measurements  reproduce  the  mesoscale  vortex  very  well.  However, 
when  the  beam  becomes  significantly  wider  than  the  core  radius, 
some  portion  of  both  the  positive  and  negative  vortex  peaks  will 
be  within  the  beam  unless  the  beam  is  entirely  to  one  side  of  the 
vortex  center.  Therefore,  one  would  not  expect  the  mean  Doppler 
velocity  value  to  maximize  until  the  beam  has  just  cleared  the 
vortex  center.  For  example,  when  the  beamwidth  is  three  times 
the  core  radius,  the  peak  is  approximately  1.5  core  radii  or  one- 
half  beamwidth  from  the  center.  Likewise,  for  beamwidth  5 and  10 
times  greater  than  the  core  radius,  the  peaks  are  at  approximately 
2.5  and  5 core  radii,  respectively.  Also,  the  wider  the  beam, 
the  greater  the  smoothing  of  the  true  tangential  velocity  profile 
and  consequently  the  smaller  the  magnitude  of  the  peak  Doppler 
velocities.  Detectability  decreases  with  range  as  the  beam 
becomes  wider  relative  to  the  vortex. 

Since  the  Doppler  velocity  profiles  appear  to  peak  at  a 
radius  of  about  one-half  beamwidth,  the  curves  in  Figure  1 were 
replotted  relative  to  beamwidth.  The  resulting  curves  (Fig.  5-4) 
vividly  protray  what  we  call  the  tornadic  vortex  signature  (TVS). 
The  peak-to-peak  diameter  is  not  significantly  affected  by  the 
size  of  the  within-beam  tornado.  However,  the  signature  ampli- 
tude, which  is  affected,  plays  an  important  role  in  TVS  detecta- 
bility. The  TVS  cannot  be  resolved  unless  the  peak-to-peak 
Doppler  velocity  shear  is  appreciably  greater  than  the  background 
cyclonic  shear  produced  by  the  parent  mesocyclone. 
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Figure  5-4.  Theoretical  change 
in  size  and  magnitude  of  toma- 
dic  vortex  signature  (TVS)  as 
radar  beamtidth  changes  relative 
to  core  radius  of  peak  tangen- 
tial velocity  in  tornado. 
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Several  other  practical  limitations  must  be  considered  when 
attempting  to  identify  tornadic  vortex  signatures  from  mean 
Doppler  velocity  measurements.  First  of  all,  when  data  are 
collected  at  discrete  azimuthal  increments,  peak  values  may  not 
be  sampled  when  the  sampling  interval  is  greater  than  one  beam- 
width.  Secondly,  Zrnic  and  Doviak  (1976)  have  shown  that  a radar 
antenna  rotating  rapidly  relative  to  the  sampling  time  has  an 
effectively  broadened  beamwidth.  Thus,  for  a given  peak  tan- 
gential velocity  and  tornado  size,  the  amplitude  of  the  TVS  will 
decrease  as  the  antenna  rotation  rate  increases. 

5.5  Tornado  Vortex  Signature  Statistics 

Objective  criteria  for  defining  a tornadic  vortex  signature 
have  not  yet  been  established.  We  are  in  the  process  of  analyz- 
ing all  signatures  and  possible  signatures  found  to  data.  A more 
thorough  search  of  the  NSSL  Doppler  radar  archives  also  is  under- 
ray. We  do,  however,  have  some  tentative  guidelines  for  defining 
a TVS.  These  guidelines  are: 

a.  An  azimuthal  shear  of  at  least  15  to  20  m s-1 
over  an  azimuthal  distance  of  1 deg  or  less  (1.25 
beamwidths  or  less). 

b.  Signature  cyclonic  (TVS  translation  removed) 
with  peak  Doppler  velocity  values  of  opposite 
sign.  (A  rare  anti  cyclonic  TVS  would  not  be 
excluded,  however.) 

c.  Anomalous  shear  region  not  more  than  about  1 km 
range  extent  (otherwise  it  would  indicate  a shear 
line  rather  than  a vortex). 

d.  Shear  region  at  least  a few  kilometers  in  vertical 
extent. 

e.  Persistent  anomalous  shear  region  at  the  same 
general  heights  for  about  10  min  or  more. 

Using  these  tentative  criteria,  we  thus  far  have  identified 
nine  signatures  in  the  1973-1975  data  set.  Since  TVS  detection 
is  a function  of  tornado  size  and  strength  as  well  as  radar 
sampling  density,  not  all  tornadoes  in  existence  during  Doppler 
radar  data  collection  have  produced  noticeable  signatures. 

The  nine  identified  signatures  are  listed  in  Table  5-2. 

Seven  signatures  had  tornadoes  or  funnel  clouds  associated  with 
them.  These  occurrences  were  either  documented  from  Storm  Data 
(Dept,  of  Commerce,  1973-1975)  or  from  damage  surveys  conducted 
by  NSSL  meteorologists.  The  remaining  two  signatures  had  no 
tornadoes  reported  with  them.  Unfortunately,  we  do  not  know  what 
occurred  at  these  rural  sparsely  populated  locations  because 
damage  surveys  were  not  conducted. 
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5.6  Summary 


Mesocyclones  associated  with  thunderstorms  and  tornadoes  can  be  identi- 
fied from  single  Doppler  mean  velocity  fields  by  applying  objective  criteria 
to  regions  of  cyclonic  shear. 


The  presence  of  a tornado  vortex  within  a Doppler  radar  sampling  volume 
also  results  in  a unique  signature  but  it  is  not  as  easily  identified.  Both 
signatures  however  do  provide  a means  for  identifying  these  circulations 
which  are  also  hazardous  to  aircraft  operations. 
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6.  GENERAL  SUMMARY  AND  CONCLUSIONS 


The  spectrum  broadness  of  the  radial  velocity  as  observed  by  Doppler 
radar  appears  to  have  a high  potential  as  an  indicator  of  thunderstorm  tur- 
bulence. In  the  76  turbulence  encounters  having  Uje  ^ 6.1  m s"'  (20  ft  s”l) 
the  turbulent  regions  were  associated  with  Doppler  radar  radial  velocity  "B" 
values  4 m s"’  within  2 km  72  times  for  a POD  = 95%.  In  addition  thunder- 
storm associated  mesocyclones  and  tornadoes  hazardous  to  aircraft  operation 
can  be  identified  using  single  Doppler  velocity  fields  by  applying  objective 
criteria. 


Utilization  of  these  parameters  in  the  air  space  system  is  dependent  on 
information  availability  and  accessibility.  Toward  this  end  the  occtant 
change  counter  method  of  estimating  mean  velocities  was  developed.  Further 
development  led  to  the  multi -moment  display  which  now  provides  on  a single 
display  in  real  time  the  reflectivity,  mean  velocity,  and  spectrum  broadness 
at  grid  locations  within  a storm.  Thus  the  goal  of  increased  safety  and 
better  utilization  of  air  space  in  thunderstorm  conditions  appears  attainable. 


TABLE  5-2 

TORNADIC  VORTEX  SIGNATURES  AND  ASSOCIATED  WEATHER  PHENOMENA 

Associated 


Date 

Location 

Phenomena 

24  May  1973 

Union  City,  OK 

Tornado,  Hail 

23  May  1974 

Yukon,  OK 

Tornado 

6 June  1974 

Tabler,  OK 

? 

8 June  1974 

Oklahoma  City,  OK 

Tornado 

8 June  1974 

Harrah,  OK 

Tornado 

13  June  1975 

Stillwater,  OK 

Tornado 

13  June  1975 

Ripley,  OK 

Funnel  Aloft 

13  June  1975 

Cushing,  OK 

Funnel  Aloft 

13  June  1975 

Kendrick,  OK 

? 

7.  RECOMMENDATIONS 
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Continued  refinement  of  the  models  and  development  for  the  use  of  Doppler 
radar  should  be  continued.  Further  development  of  the  real  time  displays  and 
automation  needs  to  be  supported.  Means  for  using  Doppler  radar  effectively 
in  the  National  Air-Space  system  should  be  developed  because  this  probe  has 
great  promise  for  mitigating  thunderstorm  associated  flight  hazards. 
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APPENDIX  A 


Doppler  Radar  Data  Processing  Program--PPI-scan 


FACTOR  - library  subroutine— controls  size 
of  all  plots  by  its  argument. 


RH  - horizontal  radius  of  influence  KM 

RZ  - vertical  radius  of  influence  KM 

X0,Y0  - SW  corner  of  grid 

Z0  - beginning  elevation  angle 

KEY  - determine  variable  to  be  fit-O-vel, 

1-std,  2-Z,  3-all 

DX,DY  - horizontal  grid  increments  KM 
DZ  - vertical  grid  increment  deg 
ITIME  - reference  time  if  horizontal  motion 
used 

DIST  - straight  line  distance  between  NRO-CIM  KM 
TAN2  - angle  NRO-CIM 


ISTN  - station  number  l-NRO(D),  2-CIM{D), 
3-NR0(57) 

ISTART,  ISTOP  - starting  and  stopping  times 
AZSTART,  AZSTOP  - starting  and  stopping  azimuth 
VBAR,  VDIR  - horizontal  motion 
MAGC  - AGC  cutoff 

loop  disregarded  if  KEY  ^ 3 


30 


I 


Setting  KNY  to  variable  to  be  fit 

0 - vel 

1 - std 

2 - Z 

3 - all  3 


Initializing  arrays 
SL  - slant  range 
V - velocity  (radial) 

STD  - standard  deviation  of  radial  velocity 
ALGZ  - reflectivity  dBZ 
X,Y,Z  - cartesian  coordinate  with  respect  to 
NRO  0°  elevation  plane 
lAGC  - power 


AGM  - Cressman  dividend 
AG  - Cressman  divisor 

MCR  - character  array  used  to  indicate  # pits 
interpolated  to  each  grid  point 


conversion  degrees  - radians 


m 


m 

m 


READ-ISTN0 

NTIME.AZM, 

ELEV.NDATE.SL, 

NGREG.V.STD, 

ALGZ,IAGC 


DDIST=0 


(D 


NTIME  - time 
AZM  - azimuth 
ELEV  - elevation 
NDATE  - Julian  date 


SL  - slant  range 
V - radial  velocity 
STD  - std.  dev. 

ALGZ  - reflectivity  dBZ 


NGREG  - gregorian  date  lAGC  - AGC  value 

SL,V,STD,ALGZ,IAGC  read  in  sets  of  32;  a cut 
has  been  made  on  STD  and  ALGZ  of  S/N  < 0. 

DDIST  « distance  from  STM  to  NRO(D) 


32 


SUBROUTINE 

GRID 


NGPU=RH/DX 


IMAX=70+NGPU 

JMAX=70+NGPU 

IMIN=1-NGPU 

JMIN=1-NGPU 


KP=1,32  WlO) 


DATA  CUTS>i%(lO) 


NGPU  - number  of  grid  points  w/i  radius  of 
influence 

i.e.  RH  = 1.5  DX  = .5 
NGPU  = 1.5/. 5 = 3 


I=ABS(X(KP) 

-X0)/DX+1.5 


J=ABS(Y(KP) 

-y0)/DY+].5 


ALPHA=ATAN 

{Z(KP)/SQRT 

(X(KP)*2+Y(KP)* 

*2)) 

i ~ 

K=(ALPHA- 

Z0)/DZ+1.5 


I > IMAX  ^ 

.OR.J  > JMAX.0R.I<1  MlK^DiS/lOj 

J^false 

I false 

I1=I-NGPU 

J1=J-NGPU 

K1=K-1 


pt(I,J,K)  - closest  grid  point  to  data  pt. 
ALPHA  - elevation  of  data  point  from  NRO(D) 
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SUBROUTINE 

PLOT 


II 

1,70 

If 

1,70 

=1,3 

AG(I,J,K)= 

AGd.J.K)/ 

AGW(I.J,K) 


false 


AG(I,J,K)= 

AG(j.j;K)/ 

lfA6W(I,J,K) 


true 

N*l,3  ' 

* 

false 

CALL  '' 
SMOOTH  ; 

true 


fals^ 


CALL 

PLOTC 


AG 


value  of  variable  calculated  for  each 
grid  point 


zWiVi 


where  W-j  = Cressman  weight 
V^-  = value  of  variable 


6th  argument  variable  from  Shuman  filter 
= .25  damps  out  2d  waves, 
i.e.  DX  = .5  1.0  km  waves  are  damped  out. 


1st  argument  contour  interval 
2nd  argument  beginning  contour 


1 


38 


41 


SO(I,J,K)=S(I,J,K) 
+.25*(S(I,J+1,K)+ 
S(I,J-l,K)-2. 
Sd.J.K) 


false 


&J=69 


true 


S(I,J,K)=SO(I,J,K) 


false 


dy)  = f(y)  + .25[f(y+d)  + f(y-d)  - 2f(y)] 

Assume  f(y)  = Ae^*^  k = ^ 

f(y)  = Ae^*^^  + .25[Ae^'‘ye^'^‘^  - 2Ae^'^^  + 
Aelkyeikd] 

= Ae^'^^Cl  + .25(e^^‘^  - 2 + e'"'^'^)] 

= Ae^'^^Cl  + .25(-2  + 2cos  kd)] 

= Ae^'^Cl  - .5(1  - cos  kd)] 

= Ae^'^^Cl  - .5(2sin2(^)] 

= Ae^'l^yCl  - sin^(M)]  k = ^ 

- «/..\n  _j_2/  kd  \ 1 I _ 1 


= ^^(y)[l  - sin^(^)]  L = wavelength 
Response  func.  = 1 - sin*(^) 

R(L)  = 1 - sinM^) 


=2.69 

=1,70 


S0(I,J,K)=S(I 

,J,K) 

+.25*(S(I+1, 

J,K) 

+S(I-1,J,K) 

S{I,J,K) 

-2.* 

Wavelength  removed  by  application  of  this 
filter. 

for  R = 0 

1 - sin"(^)  = 0 

sind^)  = 1 

M - !L  hi  ^ 

2~2*  2 * 2 ••• 

k = -p  = wave  # L = wavelength 


d TT  Stt  5ir 
L ■ 2 ’ T ’ 2 

L * 2d,  Id  , Id 


APPENDIX  B 


Data  Presentation  Procedures 


The  data  includes  Doppler  and  aircraft  data  for  1974,  1975  and  1976. 

Time  and  derived  gust  velocity  for  each  year's  penetrations  are  archived 
on  a single  tape.  Tom  Jobson  and  Phil  Bothwell  wrote  programs  to  build 
corresponding  Doppler  tapes  containing  only  that  Doppler  data  required  for 
use  in  processing. 

In  order  to  obtain  the  aircraft  position  relative  to  the  Doppler,  the 
digital  recorded  WSR-57  and/or  Doppler  reflectivity  data  were  first  plotted 
using  Calcomp  plotter.  These  plots  were  aligned  with  the  WSR-57  film  data 
in  order  to  minimize  azimuthal  or  range  variations  due  to  scope  synchro- 
drive lag.  The  aircraft  position  is  shown  on  the  film  by  the  transponder 
[also  called  Identification  Friend  or  Foe  (IFF)]  and  the  primary  return  or 
"skin-paint".  These  plane  positions  are  traced  on  the  plot.  If  a Doppler 
plot  was  not  available  for  the  same  elevation  angle  as  the  WSR-57,  the 
aircraft  positions  from  the  WSR-57  plot  were  directly  transferred  to  a 
Doppler  grid  with  no  correction. 

Using  the  set  of  points  which  represented  the  plane's  position  at 
20  second  intervals,  a least  squares  fit  of  time  vs.  range  and  time  vs. 
azimuth  is  used  to  compute  the  plane  position  at  the  midpoint  of  every  five 
second  interval  from  the  beginning  of  the  penetration  to  the  end  of  the 
penetration. 

For  each  five  second  interval  thus  determined,  the  Rough  Rider  turbu- 
lence data  are  scanned  and  the  maximum  recorded  turbulence  is  found  and 
stored. 

Next,  the  Doppler  data  (reflectivity,  radial  velocity  and  standard 
deviation)  are  read  from  the  Doppler  tape  and  stored.  From  the  Doppler 
data,  the  Laplacian  of  the  velocity  is  computed.  The  gradient  of  the  velo- 
city is  also  computed,  along  with  the  gradient  of  the  gradient  and  the 
reflectivity  gradient.  The  "radial  shear"  is  the  second  partial  derenaline 
of  the  velocity  over  three  radials.  The  "gate  shear"  is  the  second  partial 
derivative  of  the  velocity  over  three  gates  along  the  same  radial! 

For  the  aircraft  position  corresponding  to  the  mid-point  of  each  five 
second  interval,  the  Doppler  data  for  an  area  within  1 km  of  the  aircraft  is 
searched  to  find  the  maximum  reflectivity,  velocity,  spectrum  width,  velocity 
gradient,  gradient  of  the  gradient,  Laplacian,  radial  shear,  and  gate  shear. 
These  maximums  are  then  printed  out  alongside  the  derived  gust  velocity. 

For  easy  visual  representation,  a printer  plotting  program  was  developed 
which  plots  the  maximum  standard  deviation,  velocity  gradient,  gradient  of 
the  gradient,  and  the  Laplacian,  on  the  same  graph  as  the  derived  gust  velo- 
city (Y  axis)  while  each  5-second  interval  is  represented  on  the  X axis. 

For  the  average  penetration  speed  this  5-second  interval  translates  into  a 
1 km  (0.5  n mi)  interval. 


